Ferroelectric charge order stabilized by antiferromagnetism in multiferroic LuFe204. 
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Neutron diffraction measurements on multiferroic LuFe2 04 show changes in the antiferromagnetic 
(AFM) structure characterized by wavevector q = (5^5) as a function of electric field cooling 
procedures. The increase of intensity from all magnetic domains and the decrease in the 2D magnetic 
order observed below the Neel temperature are indicative of increased ferroelectric charge order. The 
AFM order changes the dynamics of the CO state, and stabilizes it. It is determined that the increase 
in electric polarization observed at the magnetic ordering temperature is due to a transition from 
paramagnetic 2D charge order to AFM 3D charge order. 



Ferroelectric polarization (FE) arises from covalent bond- 
ing between anions and cations or the orbital hybridiza- 
tion of electrons and generally excludes occurrence of fi- 
nite 3d magnetic moments and associated magnetismii 
However, competing bond and charge order in transition 
metal compounds can lead to inversion symmetry break- 
ing and multiferroic behaviour 1^ Alternatively, FE polar- 
ization may arise from frustrated charge order (CO) as re- 
ported for LuFe2 04.'^ This compound is of particular in- 
terest, as, in addition to ferroelectricity, magnetism orig- 
inates from the same Fe ions and this holds the promise 
of strong magnetoelectric coupling. The FE and mag- 
netic order (MO) takes place at and near ambient tem- 
perature, respectively, which provides the potential for 
room temperature multiferroics. These properties make 
them extremely useful in devices where the coexistence 
of both charge and spin components can be exploited and 
one property can be used to drive the other. 

LuFe204 has generated a large scientific interest since 
its multiferroic properties were reported by Ikedai^ 
LUFC2O4 adopts a FE ground state below Tco ~350 K, 
while below Tjv ~240 K the Fe magnetic moments or- 
der which enhances the FE polarization by 20%. It has 
been demonstrated that the dielectric properties can be 
changed with applied magnetic field^i^ and that pulsed 
currents can change the magnetic state^ii although the 
latter is under debate.^ These studies confirm the strong 
coupling between the magnetic and FE orders, yet the in- 
trinsic nature of the charge order suggests the charge and 
magnetic order are in principle independent 1^ For exam- 
ple, it has been recently demonstrated that the magnetic 
and charge fluctuations are largely independent in the 
frustrated magnet NiGa2S4<^ 

Various electron, x-ray and neutron diffraction (ND) 
studies have investigated the charge and magnetic or- 
ders. The crystal structure of rhombohedral LuFe2 04 
(_R3m) consists of a triangular double layer of iron ions, 
with trigonal bipyramids of five oxygen nearest neighbors 
(n.n.), in which an equal amount of Fe^+ and Fe'^"'" are 
believed to coexist at the same site (Fig[TJa)). The ob- 
servation of the (■§■-§■■§■) reflections using resonant x-ray 
Bragg diffraction (RXD) further supports the existence 
of C0,^° and we demonstrated that the Fe^"*" orbitals 
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FIG. 1. (a) Drawing^^ of the FeOs bilayer indicating small 
red Fe'^^ and large blue Fe^^ ions with magnetic moments 
oi 5 fj,B and 4 fj,B- The y/3 x \/3 unit cell is indicated in 
bold solid lines and the total magnetic moment per layer is 
indicated at the left. (b,c,d) AFM order for (b) AFE CO, (c) 
FE CO and (d) disordered CO (2D) respectively. The n.n. 
distance of iron ions is 3.65 A m a smg le layer, 3.23 A m a 
bilayer, and 6.34 A between bilayers. 



exhibit a glass like order.— 

The triangular lattice gives frustration for a 1:1 ratio of 
Fe^"*" to Fe'^'^ ions, but is stable for a 1:2 ratio where one 
species is surrounded by a honeycomb lattice of the other. 
Thus one of the layers has a majority of Fe^"*" ions and 
the other a majority of Fe^+ ions, which results in n.n. of 
the opposite species for all minority ions. Such CO gives 
rise to an electric dipole moment. While the CO state in 
the bilayers is understood, the sign of electrostatic inter- 
action between the bilayers is unclear. The spontaneous 
polarization suggests a FE ground state, where the FE 
moments of the bilayers are aligned along the c-axis. In 
contrast an antiferroelectric (AFE) ground state has been 
suggested from x-ray diffraction^^ 
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FIG. 2. ND intensity recorded in the [H H L] plane at (a) 
T=225 K and (b) T=180 K, after cooling the sample from 
T=350 K in zero applied electric field. At T=225 K a broad 
magnetic rod is observed at (ii L) characteristic of 2D mag- 
netic ordering. The peak observed at (HO) is characteristic 
of FM alignment of the FeOs bilayers. At T=180 K mag- 
netic superlattice peaks are observed at (||-§-) characteristic 
of AFM alignment of the FeOs bilayers. 



In the present study we investigate the magnetic or- 
der in LuFe2 04 with single crystal ND as a function of 
applied electric field cooling procedures. We observe a 
change in the magnetic long range order when an applied 
electric field is present at Tjy. Our results show that long 
range CO that is promoted by the electric field, is stabi- 
lized by the magnetic order. We propose that AFM order 
in the Fe bilayers facilitates the formation of long range 
CO because electron hopping between the sheets of the 
bilayer is assisted by the alignment of the Fe^+ and Fe"^"*" 
spins. 

A LuFe2 04 single crystal was prepared as reported 
previously and cut from the same ingot as the crystals 
used for the RXD study.— Pyroelectric current measure- 
ment shows a decrease of the electric polarization above 
T=220 K of ~10% and the sample becomes conductive 
due to the loss of CO above Tco ^^330 K. 

ND experiments were performed at the Wombat 
diffractometer at the OPAL facility in Sydney, using 2.95 
A monochromatic neutrons. The wavelength was cali- 
brated using AI2O3, and the efficiency of the area detec- 
tor was calibrated with a vanadium sample. The diffrac- 
tion peaks of interest were selected from the area detec- 
tor to reduce the background. We noted several other 
crystallites besides the one of this study. The selected 
single crystal was aligned with the [110] and [001] axes 
in the scattering plane and mounted between two vana- 
dium electrodes giving rise to an electric field in the scat- 
tering plane at an angle of 30 degrees to the [110] axis. 
The sample was electrically isolated since thermal con- 
tact was achieved with a sapphire rod. The electric field 
was applied in vacuum without contact to the sample but 
surface currents or shorts through a part of the sample 
cannot be excluded. The electric field was removed and 
thermal equilibrium was achieved before each ND mea- 
surement. 

We recorded ND spectra at T=225 K and at T=180 K, 
following a sequence of zero electric field cool (ZEFC) 



and electric field cool (EFC) sequences. First, the sam- 
ple was heated to T=350 K and cooled without electric 
field to obtain a ZEFC ground state. Fig. [2] shows the 
diffraction intensity along (i-lL), recorded at T=225 K 
and T=:180 K respectively. The rod of intensity along 
(iiL) observed at T=225 K is characteristic of 2D mag- 
netic order in the ab planes. Similar behavior is reported 
by Ref^. The g-dependence of the intensity is well de- 
scribed by 2D magnetic order with magnetic moments 
along the c-axis. Because of the low q range of our ex- 
periment we are not sensitive to the CO. 

At T=180 K superlattice peaks are observed at half 
integer values of L, indicating antiferromagnetic (AFM) 
order, while superlattice peaks at integer L are largely ab- 
sent. The refiection observed at (^-jO) does not change 
as function of temperature and electric field. It is be- 
lieved to be distinct from the AFM structure and will 
be discussed later. First we focus on the half integer su- 
perlattice. In contrast with ref J^, where the half integer 
peaks are attributed to ferrimagnetic alignment of Fe^+ 
and Fe^"*" moments, we observe all half integer L without 
the presence of all integer L refiections. This suggests 
that AFM and ferromagnetic (FM) arrangements of the 
bilayers occur in distinct regions of the crystal and are 
possibly formed depending on the stoichiometry of the 
sample. 

Symmetry allows for three CO domains that are char- 
acterized by propagation vectors p^=(i, i, |), Pb = {-^: 
i, |) and pc={^, |)^ and three magnetic domains, 
characterized by mi= (i, i, |), m2= (^, |, |) and 
m3=(5, "^i §)■ Fig. H^a)) shows a bilayer with p^i and 
m2. Order with and has the same motif, while 
order with p^ and mi is different. Distinct magnetic 
intensities are calculated from the structure factor, and 
the similar intensity suggests a large amount of magnetic 
and/or lattice stacking faults along the c-axis (Fig. 

Further, the sample was heated to T=350 K and cooled 
to T=225 K in an applied electric field of 350 kV/m. 
The neutron intensity recorded at T=225 K is identical 
to that of the ZEFC case. We emphasize that the 2D 
rod observed at T=225 K is magnetic while the CO is 
predominantly 3D, as observed with RXDi^i 

The electric field was removed before the ND was 
recorded, and before the sample was further cooled to 
T=180 K. Fig [3] compares the ND intensity recorded at 
T=:180 K with the ND intensity recorded after ZEFC. 
The intensity of domain D2 has increased significantly 
while the intensity of domains Dl and D3 has decreased, 
indicating a change in magnetic domain population after 
EFC to T=225 K. The width of the superlattice peaks 
along [110] is limited by the resolution of the instrument 
while it is larger than the resolution along [001] and well 
described with Gaussian peak profiles. The magnetic cor- 
relation length is taken as the inverse of the half width 
at half maximum. The fitted integrated intensities and 
correlation lengths are summarized for the three domains 
in the inset of Fig. [3] and in Table HI 

The correlation length of the majority domain after 
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FIG. 3. Neutron intensity recorded along the L) direc- 
tion at T=180 K as function of EFC procedure. The mag- 
netic state obtained after EFC to T=225 K is compared to 
the magnetic state obtained after ZEFC (bottom) and to the 
magnetic state obtained after EFC to T=180 K (top, shifted 
for clarity). The black dotted line indicates 2D magnetic or- 
der in the ah plane with moments along the c axis. The inset 
shows the integrated intensity, normalized to the total inten- 
sity at ZEFC, as function of EFC sequence. Dl, D2 and D3 
correspond to mi, m2 and ms respectively. 



EFC to T=225 K (D2) has increased, while the correla- 
tion length of the other two domains (Dl and D3) has 
decreased. This can be understood as a preferential pop- 
ulation of one type of CO domain due to the electric 
field present at Tco leading to larger correlation length 
of these domains along the c-axis. 

Finally the sample was heated to T=225 K and cooled 
in an electric field (EFC) of 350 kV/m to T=180 K. The 
electric field was removed before the ND was recorded 
(Fig. [31). Remarkably the intensity of all magnetic do- 
mains has increased significantly due to the electric field 
present at the magnetic ordering temperature. At the 
same time the intensity of the diffuse rod, observed as 
the high background in Fig. [31 has decreased. 

The correlation lengths after EFC to T= 180 K are 
largely the same as after EFC to T=225 K (TableHl). The 
additional magnetic intensity originates from magnetic 
domains of the same size, suggesting that the electric field 
assists in the nucleation of additional magnetic long range 
order. Alternatively, the additional magnetic intensity 
can originate from an increase in the magnetic structure 
factor, as discussed next. 

As now well understood, the CO in a single FCO5 
bipyramid layer consists of a 1:2 (or 2:1) ratio of Fc'^^ 
and Fe'^"'" ions. Such an arrangement is most stable for 
the Ising model on the triangular lattice. Because the 



TABLE I. The correlation lengths (in A) of the three magnetic 
domains along the c-axis recorded at T=180 K as function of 
electric field cooling procedure (errors between brackets). 

Dl D2 D3 

ZEFC 82(2) 75(3) 84(1) 

EFC to T=225 K 70(2) 92(1) 76(1) 
EFC to T=180 K 85(1) 95(1) 75(1) 



magnetic moments exhibit a strong magnetic anisotropy 
along the c axis as well as an AFM n.n. interactior^i^, a 
similar arrangement for the magnetic spins is expected, 
and indeed observed." If the CO and MO are in-phase 
(in a single layer) all Fe^+ moments are parallel and form 
one sublattice, and all Fe'^+ are antiparallel to the Fe^"*" 
moments, forming a second sublattice, see Fig. [ija) bot- 
tom layer. However, the CO and MO can be out-of-phase 
which results in a more complex structure as indicated 
in the top layer of Fig. [Ija). Indeed it has been recently 
observed from magnetic x-ray dichroism studies'^='l^ that 
the majority Fe^'^ layer exhibits CO and MO in-phase 
while the minority Fe^"*" layer has CO and MO that are 
out-of-phase. This results in a bilayer where all the Fe^"*" 
moments are aligned parallel, while the Fe'^'^ moment are 
either parallel or antiparallel (FiglH). This magnetic ar- 
rangement in the bilayer has also been reported from ND 
studies by Ref.-^. 

The spin only magnetic moments of the Fe^"*" and Fe'^''" 
ions are 4 and 5 /is , respectively. Because these mag- 
netic moments are distinct, the total moment of the mi- 
nority and majority layer is different even if the magnetic 
sequence of up and down moments is the same. The to- 
tal moment of a single layer is also different for MO that 
is in-phase or out-of-phase with the CO. Taking the spin 
only values and the magnetic sequence mentioned earlier, 
the above mentioned bilayer has = 3 /is per three 
iron ions in a majority Fe^"*" layer and A/+ — A per 
three iron ions in a minority Fe^''" layer. 

According to Hund's rules, Fe'^+ has five 3d electrons 
with spin up, and Fe^"*" has five 3d electrons with spin 
up and one with spin down. This latter 3d electron can 
jump from Fe^"*" to the Fe"^^ much easier if the spin down 
states of the latter are empty. This is the case for parallel 
alignment of the Fe^+ and Fc'^+ moments. Thus electron 
hopping in the majority Fe^"*" layer is arrested because 
the Fe^+ and Fe'^^ ions have opposite moments, while 
3d electron hopping is promoted in the minority Fe^+ 
layer because some neighboring Fe^+ and Fe'^+ ions have 
parallel moments. We also note that electron hopping 
between Fe^"*" and Fe'^"'" ions from one to the other layer 
is promoted, as their moments are parallel. This makes 
a convenient mechanism to switch between FE and AFE 
arrangements of the bilayers. As such the magnetism 
may assist the FE order to align with the electric field 
and enhance the electric polarization. 

We now turn to the interactions between the bilayers. 
Discussions are ongoing whether AFE or FE alignment of 
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the bilayers is the ground state, suggesting these configu- 
rations are close in energy. AFM bilayer interactions have 
been reported by Ref.^® while FM bilayer interaction was 
reported more recently. In the sample of interest in the 
present study the interaction is AFM, as witnessed by 
the half integer L superlattice. 

Fig [T]Jb,c) shows the AFM configuration for FE and 
AFE alignment of the bilayers. Because the total mo- 
ment of the individual layers of the bilayer are M~ — 3 
fiB and M+ = 4 /z^ per three Fe ions respectively, the 
magnetic structures of FiglUb) and[TIc) are distinct. In 
the case that the MO is 3D while the CO is 2D, the or- 
der of M~ and M+ layers is random as drawn in Fig. 
[Ijd). Such a magnetic arrangement gives rise to Bragg 
diffraction from layers with magnetic moments of 3 /is 
while 1 fiB adds to the diffuse 2D rod. We estimate that 
a charge disordered state (Fig. [TJd)) gives rise to ^^30% 
less intensity at the L/2 reflections, compared to AFM 
order in a FE or AFE arrangement (Fig. [ijb) and (c)). 

Therefore, the observed increase in intensity at all L/2 
reflections, combined with the decrease of intensity of 
the 2D rod, suggests that the electric field promotes the 
FE CO state over a 2D CO state. Furthermore, this FE 
state seems stabilized in the presence of AFM order as 
the measurements are performed after the electric field is 
removed. 

This interpretation is corroborated by the unchanged 
magnetic correlation length observed with or without 
EEC from 225 K to T=180 K. Once the CO bflayers are 
aligned by the electric field, the partial magnetic disor- 
der is removed and the AFM order parameter increases 
while the domain size remains the same. We note the 
correlation length for Dl (Table HI increased after EEC 
to T =180 K compared to EFC to T= 225 K. This possi- 
bly indicates that, in addition to the mechanism that we 
discuss, some regions may exhibit long range magnetic 
order after EFC only. 

The intensity of the diffuse rod after EFC to T=180 K 
remains significant indicating that a large part of the 
sample shows 2D magnetic order. Given the diversity 
in magnetic responses reported for LuFe204 samples it 
is not unreasonable to suggest that FM and AFM re- 
gions can coexist due to near degenerate magnetic ground 
states. We observe a strong refiection at (^^O) and a 
weaker one at (-11-3), both present at T=180 K and 
T=225 K. In addition a broad peak appears at {^^-2) 
at T=225 K. These intensities are characteristic of FM 
coupling of the bilayers and they are independent of tem- 
perature and electric field cooling procedure. 

It is tempting to attribute the increase in polarization 
observed at the magnetic ordering temperature due to 
the stabilization of long range CO by AFM regions. For a 
pyroelectric measurement the sample is cooled with elec- 



tric field present at the magnetic ordering temperature 
T TV ■ This suggests that the degree of increase of electric 
polarization is related to the amount of 2D CO present 
at Tjv, which is well below Tco- 

The temperature dependence of the hopping frequency 
reported by Xu et shows charge fluctuations well be- 
low Tco- Nagano et al.^'^ calculate that spin order stabi- 
lizes the CO through frustration of the spin system which 
results in an entropy gain, and enhances the electric po- 
larization below Tat. The present study suggests that 
the CO is stabilized with a well ordered spin system. 

Wen et a& observe that the long range CO decreases 
when the sample is cooled in applied magnetic fleld of 7 T 
along [110], suggesting that there is correlation between 
the orientation of the magnetic moments and the forma- 
tion of long range order, even in the paramagnetic state. 
This can be explained with the same phenomenon, paral- 
lel magnetic moments on neighboring Fe^"*" and Fe^+ ions 
reduces the energy needed for the electron to hop from 
the Fe^"*" to the Fe'^"'" ion. We conclude that the align- 
ment of the magnetic moments facilitates the hopping of 
electrons, promoting 3D FE order under applied electric 
field and promoting 2D CO without the presence of an 
electric field. 

We note that Wen et aB- also reports that in the case 
of reduced long range CO, below Tjv, the magnetic in- 
tensities at half-integer L are enhanced while those of 
integer L values have decreased. This is consistent with 
our observations and suggests the FM interactions be- 
tween the bilayers occur in regions with 3D CO, while 
AFM interactions occur in regions with 2D CO. 

In conclusion, we have demonstrated that the AFM 
order in LuFe204 can be changed with electric field cool- 
ing procedures. We propose that the increase of the elec- 
tric polarization at the magnetic ordering temperature 
is due to regions that exhibit a transition from 2D CO 
with paramagnetic moments above Tn to 3D CO with 
AFM moments below Tjv. While glassy orbitalsii and 
random magnetic moments give rise to a distribution of 
hopping integrals adding frustration to the charge order, 
the AFM order removes part of this frustration, causing a 
strong interaction between magnetic and charge degrees 
of freedom in this system. These findings add to the un- 
derstanding of the interaction between spin and charge, 
relevant for unconventional superconductivity, magneto- 
resistance, metal insulator transitions and other exotic 
electronic properties that arise in frustrated spin and 
charge systems. 
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